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Abstract 

A dynamical modelQ, 0, 0] of electroweak pion production reactions in the A(1232) region has 
been extended to include the neutral current contributions for examining the local Quark-Hadron 
Duality in neutrino-induced reactions and for investigating how the axial N-A form factor can 
be determined by the parity violating asymmetry of N(e, e') reactions. We first show that the 
recent data of (e, e') structure functions Fx and F2, which exhibit the Quark-Hadron Duality, are 
in good agreement with our predictions. For possible future experimental tests, we then predict 
that the structure functions Fx, F%, and F3 for (u,e) and {y,v') processes also show the similar 
Quark-Hadron Duality. The spin dependent structure functions gx and 52 of (e, e') have also been 
calculated from our model. It is found that the local Quark-Hadron Duality is not seen in the 
calculated gx and g2, while our results for gx and some polarization observables associated with 
the exclusive p(e, b'tt) and p(e, c'tt) reactions are in reasonably good agreement with the recent 
data. In the study of parity violating asymmetry A of N(e,e') reactions, the relative importance 
between the nonresonant mechanisms and the A excitation is investigated by taking into account 
the unitarity condition. Predictions are made for using the data of A to test the axial A^-A form 
factors determined previously in the studies of N(vu, (jTir) reactions. The predicted asymmetry 
A are also compared with the Parton Model predictions for future experimental investigations of 
Quark-Hadron Duality. 

PACS numbers: 12.15.Ji, 13.60.Le, 25.30.-c 



1 



I. INTRODUCTION 



In recent years, a dynamical model had been developed 0, HI to investigate electromag- 
netic pion production reactions in the region near the A(1232) resonance. The model was 
subsequently extended^ to also investigate neutrino-induced reactions (u^, pTir). Fairly 
consistent descriptions of all of the available data in the A region have been obtained. In 
this work, we further extend this model (called the SL model ) to address two questions of 
current interest : (1) Will Hadron-Quark Duality first observed by Bloom and Gilmanj^] in 
inclusive p(e, e') be also seen in the neutrino-induced N(u, e) and N(u, v') reactions ? (2) 
How can parity violating asymmetry of inclusive N(e, e) reaction be used to improve our 
knowledge about the axial N-A form factor ? Experimental data for testing our predictions 
can be obtained at new neutrino facilities^ and Jefferson Laboratory (JLab)0,0|. 

Our first step is to construct neutral currents within the SL model. This then allows us to 
extend our previous calculations [3] of N(u^, /i~vr) reactions to also predict N(u, u'tt) reactions 
for examining some questions concerning the Quark-Hadron Duality. Simply speaking, an 
inclusive lepton scattering observable exhibits Quark-Hadron Duality when an average of 
this quantity over an appropriately chosen scaling variable in the resonance region is close 
to that in the Deeply Inelastic Scattering (DIS) region. More details on this subject can be 
found in a recent review by Melnitchouk, Ent and KeppeljH). The Quark-Hadron Duality 
was first observed^] for the p(e, e') structure function F 2 = ujW 2 . Recent experimental 
data@, 03 for both F 2 and Fi = m^Wi have further confirmed more quantitatively this 
interesting observation. Theoretical attempts in understanding the Quark-Hadron Duality 
within QCD were alre ady made 11 . 12L 13l | in 1970's. More recent works are reported in, 



for example, Refs. 0, H H |31iJllJwithiii the Standard modelQ, it is natural to 



ask whether the Quark-Hadron Duality should also exist in the neutrino-induced N(u, e) 
and N(u, u') reactions (from nowon, v and e denote the neutrinos and charged leptons, 
respectively, of any generation within the Standard Model). In the absence of necessary 
neutrino data, this can be explored theoretically in the A region using the extended SL 
model presented in this paper. Obviously, we can only explore the 'local' Quark-Hadron 
Duality in the A region. 

We will also report on our investigations of the spin dependent structure functions g\ and 
g 2 of (e, e') processes. To see the accuracy of our model in predicting these quantities, we 
first show that some recent data0, 2^] of the polarization observables associated with the 
exclusive p(e, e'n) and p(e, eV) reactions agree well with the predictions from the considered 
SL model. We then find that g\ calculated from our model agree reasonably well with 
the recent data[23|, but do not show Quark-Hadron Duality when compared with the DIS 
data 24l . The calculated g 2 is also found to be rather different from the DIS data which can 
be described reasonably well by the Wandzura-Wilczek formula (25^. 

Our second task in this work is to apply the extended SL model to address some questions 
concerning the N-A form factors. It has been well recognized that these form factors are 
important information for testing current hadron models and also lattice QCD calculations 
in the near future. The vector parts of the N-A form factors in the Q 2 < 4 (GeV/c) 2 
region have been rather well determined by analyzing very extensive high precision data of 
electromagnetic pion production reactions, as reviewed in Ref. |26| . On the other hand, the 
axial vector N-A form factor is not well determined mainly because of the short of high 
precision data of neutrino-induced pion production reactions in the A region. For example, 
the axial N-A from factor determined in Ref. [3J] is rather different from the one determined 
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previously [13]. This leads to uncertainties in interpreting the axial N-A form factor in terms 
of the hadron structure calculations such as those reported in Refs. (iil Ejj 13(1 l3l| . The 
situation will be improved when the data from new neutrino facilities become available in 
the near future. Alternatively, progress can be made by following Refs. 
H, 3^, to investigate the parity violating asymmetry of inclusive N(e, e') reactions. This 
polarization observable is due to the interference between the electromagnetic currents and 
neutral currents and hence can be used to explore the axial N-A form factor. Experimental 
effort in this direction is being made 6] at JLab. To facilitate this study, we have applied the 
extended SL model to explore the dependence of the parity violating asymmetry of N(e, e') 
on the axial iV-A form factor. 

In section II, we will specify the considered models of electromagnetic currents (em), 
weak charged currents (cc), and weak neutral currents (nc). We first recall their forms in 
the Standard model [20] and then specify how such currents are defined in terms of hadronic 
degrees of freedom of the SL model. In section III, we give expressions of the inclusive cross 
sections and indicate how the structure functions are calculated within our hadronic model 
and Parton Model. Section IV is devoted to present results for studying the Quark-Hadron 
Duality in the A region. The predicted N(e, e') structure functions will also be compared 
with the recent data. The results for the parity violating asymmetry of the inclusive N(e, e') 
will be given in section V. A summary is given in section VI. 



II. MODELS OF ELECTROWEAK CURRENTS 

We first recall the electroweak currents defined in the Standard Model[2(j. In the con- 
sidered A excitation region, we can eliminate heavy W and Z bosons and keep interactions 
involving only up (u) and down (d) quarks. The interaction Lagrangian for our study can 
then be written as 

L eff = -Vbmie-fe + j^Ap 

-^7m(1 - 7b> + e(2^ 7 M - 2g e Al ^)e]f nc , (1) 

where a = 1/137, G F = 1.1664 x 1(T 5 GeV~ 2 , = -1/2 + 2 sin 2 6 W , g\ = -1/2, is the 
photon field, e and v are the field operators of the charged leptons and neutrinos, respectively. 
The Weinberg angle Qw is known empirically to be sin 2 9w = 0.231 and V ua - = 0.974 is the 
the Cabibbo-Kobayashi-Maskawa (CKM) coefficient. The electromagnetic current (j^ m ), 
weak charge current (J£.) and weak neutral current (j^ c ) carried by u and d quarks can be 
written as 

j? m = lu^fu-ld-fd, (2) 

f cc = firfil-^d, (3) 

JL = \[*flX - 7b)« - dY(l - l 5 )d] - 2 sin 2 9 w f em . (4) 

With the simplification that only u and d quarks are kept, it is well known (37l| that the 
above currents can be classified according 'strong' isospin. We thus can write 

Jem 3 ' v isoscalar ' \ / 
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fee = (Vr + iV?) - {At + iA$) , (6) 
3& = (Vf - ^) ~ (A^ - %A£) , (7) 
f nc = (1 - 2 sin 2 6 w )f em - Vf soscoXar - A% , (8) 

where the isospin components of the vector (V) and axial vector (A) currents are defined as 

Vr = tfZq, (9) 

Af = rfft^q, (io) 

Vtsoscalar = ^Vff ■ (H) 

Here we have defined an isospin doublet field operator q = g + 7° = (u + , d + ) / y° and 7$ is the 
usual Pauli operator. 

The above expressions will allow us to calculate electroweak structure functions of deeply 
inelastic (e, e'VfV, e) and [y, u') processes within the Parton Model, as explained in, for 



example, Ref. [41|. However, they can not be used directly for investigating meson produc- 
tion reactions in the resonance region where the perturbative QCD is not applicable. At 
the present time, the most tractable ways for investigating these reactions are in terms of 
hadronic degrees of freedom. The starting point is a hadronic effective Lagrangian con- 
strained by the symmetry properties of the Standard Model. Accordingly, the resulting 
electroweak currents have the forms of Eqs.(5)-(8), but are written in terms of hadronic 
field operators. This can be achieved by using the standard effective chiral lagrangian 
met hods [5(1 H2I 43[. In the SL model developed in Refs.[l|, 0, 0|, the constructed elec- 
troweak currents are expressed in terms of the field operators of the nucleon (N), Delta (A), 
pion (7?), omega meson (u), and rho meson (p). They can be written as 

V» ■ t?„ = A%% - ^—a^d^} ■ T -N + 2±N>f lf $ ■ t}N x 7? 

+ [7? x cPt?] ■ v„ - 9 -^Le a ^{d a v»] ■ n[d»iu 5 ] - xA^f ■ v v T^N , (12) 



m 



T 



A^ ■ Vfj, = g A N^ l5 - • v^N - Fd»7? ■ v„ - f p7TA (^ x vf) ■ ^ + A U T ■ v^N , (13) 

where is an arbitrary isovector function, F = 93 MeV is the pion decay constant, and 
qa = 1-26 is the axial coupling strength of the nucleon. The other parameters as well as the 
vertex T^ v and isospin operator T of the iV-A transition are given in Refs.jll, H 0. With 
the electromagnetic current j£ m given explicitly in Ref.jjJ, the relation Eq.(5) and Eq.(12) 
can be used to get V^ oscalar and hence j£ c of Eq.(8) is also completely determined within the 
SL model. The neutral current j% c is needed to extend the calculations of Refs. 0, 0, HI to 
also study [y, z/'tt) reactions and the parity violating asymmetry of the (e, e') reaction. 

The electroweak currents defined by Eqs.(5)-(8) and (12)-(13) must be supplemented by 
additional effective Lagrangians describing the haronic interactions, such as nNN, ttNA, 
pNN, pTTTr, and uNN couplings, to calculate the meson production amplitudes within the 
SL model. The details of such a dynamical approach have been given in Refs. [l], 0, HI . Here 
it is sufficient to just illustrate schematically the basic meson production mechanisms of the 
SL model. In Fig.l, we show the constructed non-resonant pion production mechanisms. 
For the em current contributions, the waved-line is photon and all diagrams contribute. For 



4 



cc and nc current contributions, the wave-line is the vector current V or axial vector current 
A. Like the em case, all terms in Fig.l contribute to nc amplitude but with different isospin 
weighting factors defined by Eqs.(5) and (8). For the cc contributions, all terms contribute 
to the amplitudes induced by the vector current V except that the p exchange in (c) should 
be excluded because it is an isoscalar interaction. The cc's amplitude due to axial vector 
current A contains only mechanisms (a), (b), (e) with only p exchange, and (f). In both cc 
and nc cases, the pion pole term due to the second term —Fd^Tr-v^ of axial current Eq.(13) 
has additional contributions, as explained in Ref . [j| • 

The A excitation mechanism is illustrated in Fig.2. For the later discussions in this 
paper, here we recall that the matrix element of the axial iV-A coupling IT in Eq.(13) can 
be written as 

< Pa | IT I p N > = d 1 (q 2 )g^ + ^P-Pa^g^ - g^q u ) 



m 



N 



2\ 



,2 



WN-i-^^PafWB, (14) 



my m N 

where p& and Pn are the momenta of the A and N, respectively, q = Pa — Pn, and P = 
(pa +Pn)- The strengths of the form factors at Q 2 = —q 2 = are identified with the quark 
model predictions 28[ and are found to be di(0) = P > i^^9a with = 1.26, D\ = 0.969, 
D 2 = —0.180, D 3 = 0.836, and c? 4 = 0. In Ref. 0], it was found that all of the available data 
of N(u/j,, p~7f) data in the A region can be well described if the g 2 -dependence of each form 
factor in Eq.(14) is taken to be (Q 2 = —q 2 ) 



di, 2 (Q 2 ) = d h2 (0)Gi A (Q 2 ) 

) m N 
Q2 _|_ m 2 



d 3 (Q 2 ) = d 3 (0) ™ N , G A N AQ 2 



with 

Gi A (Q 2 ) = (1 + aQ 2 )exp(~bQ 2 )G A (Q 2 ) , (15) 

where a = 0.154 (GeV/c)~ 2 , b = 0.166 (GeV/c) 2 , and G A (Q 2 ) = 1/(1 + Q 2 /M 2 A ) 2 with 
Ma = 1-02 GeV is the nucleon axial form factor |44|. The corresponding vector form factor 
G\ N (Q 2 ), which is associated with Yy' of Eq.(12) and has been well determined |2j by 
analyzing the data of pion electroproduction, has the same form as Eq.(15) except that 
Ga(Q 2 ) is replaced by the usual proton form factor G P (Q 2 ) = 1/(1 +Q 2 /A 2 ) 2 with A 2 = 0.71 
(GeV/c) 2 . 

We further mention that the electroweak pion production amplitudes (denoted as J from 
nowon) calculated within the SL model can be written schematically in the center of mass 
(cm.) frame of the final irN system as 

(n . r .) Yj N (k,W)Y A>aN (q,W) 
J,N,aN{k, q, W) = J v nN J N {k, q, W) + _ — _ — — , (16) 

where a = 7, W or Z bosons, q and k are the relevant initial and final momenta, ^(W 7 ) is 
the A self-energy. The non- resonant amplitudes J^ N r J N (k,q,W) are calculated from 

J-KN,'aN(k, 9> W) — v wN,aN + t 7r N,TrNG 7T N(W)v w N,aN , (17) 
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where v % N, a N are the non-resonant amplitudes illustrated in Fig.l, t^N^N is the non-resonant 
tcN scattering amplitude, and G n N(W) is the 11N propagator. Note that the second term 
of Eq.(17) is the consequence of the unitarity condition. This is neglected in the previous 
investigations |35L liR I39L l4o| of the effects due to the non- resonant amplitudes on the par- 
ity violating asymmetry of N(e, e') reaction. Furthermore, the non-resonant mechanisms 
considered in those works also differ from what are illustrated in Fig.l and well tested in 
extensive investigations of exclusive N(j,7r), iV(e, eV), and N iy„ . jU.~7r) reactions. For ex- 
ample, the vector meson exchanges are not included in Refs. |38L l39j | and pseudo-scalar irNN 
coupling, instead of the pseudo- vector coupling used in the SL model, is used in Ref.pZj]. 

The resonant term in Eq.(16) is defined by the dressed vertices which contain the influence 
of the non-resonant interactions as given by the following equation 

f aM<i, w ) = FaMq) + J dkr A M k ) G Mk, w)j ( ^ N %(k, q , w) , (is) 

where r a,on(q) is the bare N-A form factor. The second term in Eq.(18) is commonly called 
the meson cloud contributions to the N-A transition. Explicit calculations of these meson 
cloud effects also mark an important difference between this work and all of the previous 
investigations [i^, HE 13 HE IIE |3?], HE HE H(j of parity violating asymmetry of N(e, e') 
reactions. We emphasize that Eqs.(16)-(18) satisfy the unitarity condition which is essential 
in interpreting the meson production data, as explained in Refs. [3, 0, H| as well as in many 
works reviewed in Ref. \2t 



III. CALCULATIONS OF STRUCTURE FUNCTIONS 

With the Lagrangian Eq.(l), the formula for calculating the exclusive cross sections for 
(e, e'n) and (v, en) processes are given in details in Refs.Jl 0, ■ The formula for calculating 
(z/, u'n) cross sections can be easily obtained from Ref. [3] with minor modifications. In this 
paper, we focus on inclusive processes. Their cross sections can be more simply written in 
terms of structure functions. 



Following Ref. [4 11] . the symmetry properties require that the lepton scattering structure 



functions W\, W2 and W3 are in general related to the hadron tensor by 

W i» = - Wig »» + I^-pV + i^- e ^p a qp , (19) 

where the convention e 0123 = —1 is chosen. The structure functions are functions of two 
independent invariant variables. One usually chooses Q 2 and the invariant mass W = \fs = 

(p + q) 2 in the resonance region, but chooses Bjorken scaling variable x = Q 2 /(2p ■ q) and 



Q 2 in the deeply inelastic region. Here q and p are the momenta of the exchanged boson 
(7, W or Z) and the initial nucleon, respectively. The unpolarized inclusive N(e, e') cross 
sections can then be written as 

da Aa 2 E' 2 r TTr 2 9 Tir 2 9, 

with 

Q 2 = -q 2 = AEE' sin 2 ^, (21) 
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where E and E 1 are the incident and outgoing lepton energies, respectively, 9 is the lep- 
ton scattering angle with respect to the incident lepton . For inclusive [y, e~) and [p, e + ) 
processes, we have 

- ' 1 £' 2 [2W\sin 2 ~ + W/ 2 cos 2 -±W 3 ^^sin 2 ~], (22) 



dE'dto' 2tt 2 1 2 2 m N 2 

where the sign in front of W3 is + (-) for (v, e~) ((z/, e + )). For {y, is'), the cross section 
formula is the same as Eq.(22) except that the factor | V U( i | 2 is removed. 

Within the hadronic models, such as the SL model considered in this paper, the hadron 
tensor W^ u in the region near the A excitation can be calculated from summing all contri- 
butions from a(q) + N(p) — > 7r(k) + N(p) with a = 7, W, Z. We can write in general the 
hadron tensor W^ v for electroweak pion reactions as 

= J2<f\jz\i>*<f\r a \i>, (23) 

where a = em, cc, nc denotes the considered current, and we have introduced concise nota- 
tions 

E = E EM 3 ^ + q- P '-k)^, (24) 



f rn N 



and 



< / I ^ I i > 



m N 



(2tt) 3 \ E N E' N 2E m 



<k(t 7T ),p'(s',t')\JZ(q)\p(s,t)>, (25) 



Here (s, t) denote the z-components of the nucleon spin-isospin, t w the z-component of the 
pion isospin, and < k(t n ),p'(s' , f)\J£(q)\p(s, t) > is the a+N — > tt+N amplitude with a = 7, 
W, Z for a = em, cc, and nc, respectively. In the dynamical approach of Refs. 0, 0, , the 
current matrix element Eq.(25) has the form of Eq.(16), consisting of a non-resonant term 
and a resonant term. 

From Eqs.(19) and (23), one can derive the expressions for calculating the structure 
functions W\, W2, and W3 from the current matrix elements defined by Eq.(25). It is 
convenient to calculate these structure functions in the center of mass (cm.) frame of the 
initial aN and final irN systems. The direction of the momentum-transfer is chosen to 
be the quantization z-direction : i.e. q = (u c , 0, 0, | q c |) for the exchanged boson and 
p = (E^, 0,0, — I q c |) for the initial nucleon. We then have for the contribution from 
current j a with a = em, cc, nc, 

wi = ~E[i</i j aM>r+i</m*>r], (26) 



2 

^2 = ^J2il(\<f\Ja\i>\ 2 +\<f\J y a\i>\ 2 ) + ^\<f\Ja\i>\ 2 ), (27) 



^3 = -^jEM< / I J a I i >< / I J l I < >*] » ( 28 ) 

where J a = J° + ^J a • Q, and q is three- momentum transfer in the laboratory frame (i.e. 
q = (co>,q)). In practice, the calculation of any term of the structure functions defined by 
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Eqs.(26)-(28) can be obtained from appropriate combinations of the following integrations 



C, P = £[< / I Z I i >< f I Jp I i > 
1 /" ,„ \k c \m N 



2 E / ^i^f < *(OpV,0I^(*)Ip(M) >< fc(^V,OTOb(M) >* • 

(29) 



S .t' .t-JT 



were a,/3 = em, cc, nc, and k c is final pion momentum in the final irN cm. system. 

We will also examine the spin dependent structure functions of (e, e'). They are 
defined |4fl| by writing the hadron tensor for a polarized target with spin vector (S 2 = 
—l,p • S = 0) as 



where 



W"" = Wg" + , (30) 



W? = wA-r + q ^f) + ^ - v -f^w - ^V) , (si) 

WT = i- q -[9iS, + gtiSp - ^Pp)) . (32) 

p • q p ■ q 



With some derivations, one can show that 



91 = yt^ [ $ xv ° + xxv ^ (33) 

up- 

92 = —^[X*-X*\, (34) 

1 + ^2 



with 



Xy° = -m^Im^ \ s=s(x) ] , (35) 
| q | 

X xy = -m N Im[lZ,em\s=s( z) ], (36) 

where I^mem \s=s(i) is the same as that defined in Eq.(29) except that the initial nucleon pro- 
jection s is not summed over but is fixed in the chosen direction defined by S(i). S(x) (S(z)) 
means that the initial nucleon spin is polarized in the direction perpendicular (parallel) to 
the incident electron direction. 

For investigating Quark-Hadron Duality, we would like to compare the structure functions 
calculated from using the above formula for hadronic models with those calculated from the 
quark currents Eqs.(2)-(4) in the deeply inelastic region. With the standard definitions 

F^x^Q 2 ) = m N Wi{x,Cf), 
F 2 (x,Q 2 ) = uW 2 (x,Q 2 ), 

F 3 (x,Q 2 ) = luW 3 (x,Q 2 ), (37) 



where x = Q 2 /(2p ■ q) = Q 2 / '(2mjva>), the Parton Model gives |4 11] (keeping only the contri- 
butions from the u and d quarks in the considered A region) 
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A. p(e,e') 



F 2 {x, Q 2 ) = x[-{u{x, Q 2 ) + u{x, Q 2 )) + -{d{x, Q 2 ) + d{x, Q 2 ))} , (38) 
F^x.Q 2 ) = F 2 (x,Q 2 )/2x, (39) 
9i(x,Q 2 ) = ^-{uHx,Q 2 )-u l (x,Q 2 )+uHx,Q 2 )-ut(x,Q 2 )] 

lo 

+^(x, Q 2 ) - d l (x, Q 2 ) + d\x, Q 2 ) - d l (x, Q 2 )} . (40) 

lo 

B. p(iy,e~) 



F 2 (x, Q 2 ) = 2xF 1 (x, Q 2 ) = 2x{d{x, Q 2 ) + u{x, Q 2 )) , (41) 
F 3 {x,Q 2 ) = 2{d{x,Q 2 )-u{x,Q 2 )). (42) 

C. p( v y) 

F 2 (x,Q 2 ) = 2xF 1 (x,Q 2 )=x[^ + (^-^sm 2 6 w ) 2 ](u(x,Q 2 )+u(x,Q 2 )) 

+x[\ + {~\ + I sin 2 e w ) 2 }(d(x, Q 2 ) + d(x, Q 2 )) , (43) 
Fz{x,Q 2 ) = [i-^sin 2 ^](n(x,Q 2 )-^(a;,g 2 )) 

sin 2 6 w }(d(x, Q 2 ) - d{x, Q 2 )) . (44) 

Here q(x, Q 2 ), q(x, Q 2 ) with q = u,d are the parton distribution functions (PDF) determined 
from fitting the data in deeply inelastic region. In Eq.(40) q^(x, Q 2 ) are the spin dependent 
parton distribution functions. The structure functions for the neutron target can be obtained 
from Eqs.(38)-(44) by interchanging the u and d parton distribution functions. In actual 
calculations, the strange quark and sea quark contributions are included, but are found 
to be very small in the considered A region. Thus our results presented below are from 
Eqs.(38)-(44). This is consistent with the considered hadronic SL model which also neglects 
any possible reaction mechanisms involving intermediate strange hadrons. 



IV. QUARK-HADRON DUALITY 

We now turn to exploring the Quark-Hadron Duality in the A excitation region. This is 
done by comparing the structure functions calculated from Eqs.(38)-(44) using the CTEQ6 
parton distribution functions [45] with those from Eqs. (26)-(28 ) and (33)- (36) using the 
hadronic SL model described in section II. It is well known|l3l. l46l |47| that more quantitative 
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tests of Quark-Hadron Duality need to include target mass corrections. Furthermore, the role 
of the higher-twist effects must be better understood. For simplicity, we will not take such a 
more involved procedure and will only compare all results from the SL model with the Parton 
Model predictions at Q 2 = 10 (GeV/c) 2 . Thus our goal here is more qualitative. We will 
focus on exploring whether the neutrino- induced reactions show the similar Quark-Hadron 
Duality observed in (e, e'). Furthermore we will also consider the spin dependent structure 
functions and parity violating asymmetry of p(e, e) to which the procedures for inclu ding 
the target mass corrections have not been developed. We follow the usual criterion|4 la. 1 1 31] 
that the local Quark-Hadron Duality is seen if the the predictions from our hadronic model 
are "oscillating" around the predictions from Parton Model such that their averaged values 
could be very close after target mass corrections are included. 

Following the previous works, as reviewed in Ref.|8j, we present the calculated structure 
functions as functions of the Nachtmann scaling variable defined by 

2x 

* = l + (l + Ax 2 m 2 N /Q 2 y/ 2 (45) 

where x = = 2m N u * s ^ e Bjorken scaling variable. One can show k| that £ = p+ /p + 
is the fraction of the plus light-front momentum of the nucleon carried by the struck quark 
in the infinite momentum frame. The use of this variable includes some of the target mass 
corrections, as discussed in Ref . [l3j| . We note here that in the £ < about 0.7 region the 
Parton Model results at Q 2 = 10 (GeV/c) 2 will correspond to W > about 2 GeV which 
is much larger than the considered A region (1.1 GeV < W < 1.4 GeV). For example, for 
Q 2 = 1.5 (GeV) 2 the A peak ( W = 1.232 GeV) occurs at £ = 0.57 which corresponds to 
W = 2.8 GeV at Q 2 = 10 (GeV/c) 2 . 

First we consider the p(e, e') processes. As reported in Refs.0, [l(|, the recent data of 
structure functions F% and F 2 from JLab have further established the Quark-Hadron Duality 
in the entire resonant region. This is illustrated in Fig. 3 along with the results calculated 
from the SL model (solid curves near £ ~ 0.6 ) in the A region. From now on, we will only 
consider the data in the A region. In Fig. 4 we compare our hadronic model calculations at 
Q 2 = 0.7, 1.5,2.5,3.5 (GeV/c) 2 (solid curves, from left to right) with these data. Each solid 
curve covers the same A region with 1.1 GeV < W < 1.4 GeV. We see that the predictions 
from the employed hadronic model (solid curves) agree well with the data and oscillate 
around the Parton Model predictions (dashed curve). 

To pursue further, it is natural to ask, from the point of view of Standard Model, whether 
the Quark-Hadron Duality observed in p(e, e') should also be seen in neutrino-induced pro- 
cesses. Experimental data for such investigations are still absent, but could be obtained at 
new neutrino facilities in the near future. To facilitate these developments, we here present 
predictions for p(u,e~) and p{v,v') processes. Since our model allows us to also predict 
the structure functions for the neutron target and hence we will also provide predictions of 
^7=0 = (Fi )P + Fi )fl )/2 for the isospin 1 = deuteron-like target. Calculations for general 
nuclear targets, such as those considered in Ref. 48||, are beyond the scope of this work. 



Our predictions for structure functions are shown in Fig.5 for (e,e'), Fig.6 for {y, e~), 
and Fig. 7 for (u, v'\ Clearly all cases exhibit the similar feature of Quark-Hadron Duality. 
Experimental confirmations of our predictions (solid curves) shown in Figs. 5-7 will be useful 
for making further progress. If they are confirmed, our next step is to examine the local 
Quark-Hadron Duality quantitatively by including the target mass corrections [12j, |4y, |47[ 
and considering the available information about the role of the higher twist effects. 
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As reviewed in Ref. |8[, an another criterion of the Quark-Hadron Duality is that the 
structure functions in the resonance region should slide along the Parton Model predictions 
as Q 2 increases. This is clearly the case in Figs. 5-7 for Q 2 < 4 (GeV/c) 2 . In Fig. 8 we further 
show within the employed SL hadronic model that this should also be the case up to rather 
high Q 2 = 20(GeV/c) 2 where the high precision data for the A region could be obtained 
from the experiments with 12-GeV upgrade of JLab. The existing SLAC data(49l. 5pl | do not 
have enough high accuracy for investigating the A at Q 2 > about 6 (GeV/c) 2 . 

We now turn to discussing the spin dependent structure functions g± and gi- We first note 
from Eqs.(33)-(36) that these functions depend on the interference terms /|^ em oc Je m Jem 
and Iemem ^ JemJem- Such terms also determine various polarization observables of the 
exclusive p(e, eV) and p(e, eV) precesses, as discussed in, for example, Ref. |5l[ . It is therefore 
important to test the SL model predictions against the recent data from such polarization 
measurements 0, 2^\. It is clear from Figs. 9-12 that the SL model can describe the JLab 



data very well and certainly can be used here to investigate the spin dependent structure 
functions. The details of these comparisons can be found in Refs. |2lL 22| and also in the 
captions of Figs. 9-12. 

Our results for g\ and #2 of p(e, e') are shown in Fig.13. In the left side of Fig. 13 we again 
see that our predictions for g\ (solid curves) agree reasonably well with the data(23j|. (Note 
that the determination of these g\ data from the polarization data of Ref. [23| involved some 
model dependent input To get the Parton Model predictions using Eq.(40), we need 
to known the spin dependent parton distribution functions. However the determination of 
such information in the considered large x region is still in developing stage and therefore no 
attempt is made here to do calculation using Eq.(40). Instead, we assume that the Parton 
model predictions can be identified with the dashed curve which is from fitting the data of 
gi measured in deeply inelastic scattering [53|. The left side of Fig.13 then indicates that 
gi clearly does not exhibit the local Quark-Hadron Duality. Such a disagreement was in 
fact expected |54T| by considering the constraints imposed by the Ellis- Jaffe integraljH^] and 
Drell-Hearn-Gerasimov sum rule 56[. 



It is known [41| that the naive Parton Model, considered here, can not predict gi- There 



are however some #2 data|57( from deeply inelastic scattering which can be described well 



by using the Wandzura-Wilczek formula [25j 

g™(x, Q 2 ) = [ f 1 ^(y, Q 2 )] - 9l {x, Q 2 ) . (46) 
Jx y 

The dashed curve in the right side of Fig.13 is obtained from the above formula using the 
dashed curve for g\ in the left as the input. Clearly it disagrees with the predictions (solid 
curves) from the considered hadronic SL model. Unlike the result for g±, there is no simple 
explanation of such a similar breakdown of local Quark-Hadron Duality. 

The Quark-Hadron Duality for the spin dependent structure functions of p(u, e) and 
p(v, v') are not investigated here since there are no corresponding PQCD calculations or 
data from deeply inelastic scattering data to compare with. 



V. PARITY VIOLATING ASYMMETRY 



The formula for calculating parity violating asymmetry A of N(e, e') reactions have been 
given in literatures 32). For our purposes, we write it in terms of structure functions defined 
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in section III. With some straightforward derivations, we have 



da(h e = +1) - da(h e = -1) 

(47) 



da(h e = +1) + da(h e = -1) 
Q 2 G F N 



y/2(Aira) D ' 



with 



N = cos 2 d -W% m ~ nc + 2 sin 2 d W{ m ~ nc + sin 2 ^(1-4 sin 2 W ) E + f Wf" , (48) 

Z Z Z 1V± 

D = cos 2 -W% m - em + 2sm 2 -Wf m - em . (49) 

Here the electromagnetic structure functions y/^ m ~ em are calculated from Eqs.(26)-(27) with 
a = em. The interference term \yem-nc between the electromagnetic and neutral currents 
can also be calculated from Eqs.(26)-(28) with the following replacements 

l< / I 4 \i >\ 2 -> < / I Jim \i >< f I 4c \i >* 

+ < / 1 r nc i i >< / 1 r em | % >* , 



< / 1 |>< / 1 jy | « >* - < / 1 j: m | i >< f | j n ^ c | i >* 

+ <f | r nc | i >< / 1 jy em | < >* . 

Obviously these quantities can also be calculated from integration 1^ defined in Eq.(29). 

We next focus on the first two terms of iV of Eq.(48) which depend on W^~ nc '. Here 
we note that W^~ nc is the symmetric part of the hadron tensor defined in Eq.(19). Hence 
they can only have the contributions from the vector parts of neutral currents because of 
the vector structure Eq.(5) of the electromagnetic current. To see this more clearly, we use 
Eq.(8) to write the neutral current as 

JL = Vnc - K , (50) 

with 

Kc = (1 " 2 sin 2 9 W ) J» m - Vf soscalar . (51) 

Obviously the [—A%] part of Eq.(50) will not contribute to W^~ nc and hence the first two 
terms of iV of Eq.(48) can then be written as 

cos 2 9 -W* m - nc + 2 sin 2 B W{ m ~ nc = (1 - 2 sin 2 ^)(cos 2 ^V 2 e "™ + 2 sin 2 B W{ m ~ em ) 
z z z z 

+ (cos 2 9 -W^ m -' s + 2 sin 2 6 -W* m ~ is ) 

= (1 - 2 sin 2 9 W )D + (cos 2 -W* m - ls + 2 sin 2 -W{ m ~ is ) , 

z z 

where D is defined in Eq.(49) and Wt m ~ %s is the same as m^ m ~ nc with J^ c replaced by the 
isoscalar vector current [— V£ oscalar ] defined in Eq.(51). Thus only the non-resonant isospin 
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1=1/2 ttN state contribute to W^ m ts . With the above relation, Eq.(47) for the asymmetry 
can then be written as 

A = — *^-[2 - Asm 2 W + Ay + A A \ 
V 2Ana 

= -Q 2 [8.99 x 10~ 5 ](1.075 + A y + A A ) , (52) 
where Q 2 is in unit of (GeV/c) 2 and 

cos 2 § Wr~ is + 2 sin 2 f Wr~ is 

A v = Y) ' ( 53 ) 

Aa = sin 2 1(1- Asm 2 6 w )^±fwr- nc (M) 

In Eq.(52) we have evaluated (2 — 4 sin 2 9yy) ~ 1.075 which is a model-independent constant 
and is the main contribution to the asy mmetr y . This term was identified by Cahn and 



Gilman[32] and later investigators |33j, |3J, 35, 36, 3jJ |38|, |3g |4 



We now note that the term A^ defined by (54) depends only on W^ m ~ nc which is anti- 
symmetric as defined in Eq.(19). Because of the vector structure of electromagnetic current 
Eq.(5), W£ m ~ nc can only have the contributions from the axial vector parts of neutral cur- 
rents. Clearly, the information about the axial vector N-A form factor, which is the quantity 
we hope to learn about as discussed in section I, is isolated in Aa through its dependence 
on W^ m ~ nc . Therefore it is important to identify the region where A^ is much larger than 
Ay such that the extracted axial N-A form factor has less model dependence. 

As mentioned above, only the non- resonant isospin I = 1/2 irN final states contributes 
to W^ m ~ ls . Thus we expect that Ay is weaker than A^ near the A resonance. This is 
illustrated in Fig. 14 for incident electron energy E = 1 GeV. We see that A^ (solid curve) 
is indeed much larger than Ay at energies near W = 1.232 GeV of the A position. This 
is also the case for other typical electron kinematics which can be conducted at JLab, as 
illustrated in Fig. 15. The contribution from A A (solid curves) are clearly much larger than 
Ay (dashed curves). Therefore the asymmetry data at resonance position W = 1.232 GeV 
can be used to extract the contributions from the axial vector currents. 

We now focus on examining how A^ depends on the axial N-A form factor. Recalling 
Eqs.(16)-(18) for the current matrix elements, A^, like other observables, contains contribu- 
tions from the resonant and non- resonant amplitudes. As seen in Fig. 16, the non-resonant 
contribution (dotted curve near the bottom) is very small compared with the full model pre- 
diction (solid curve). This agrees with previous investigations 35|, 38|, 3^, 4(J, although those 



earlier works neglected the unitary condition as discussed in section I. The large difference 
between the solid and the dotted curves in Fig. 16 indicates that A^ is a useful quantity 
for extracting the dressed axial N-A form factor defined by Eq.(18). When the pion cloud 
effects (second term of Eq.(18)) on the iV-A transition is turned off, we obtain the dash 
curves. This is consistent with our previous findings in Refs. [3, 0, HI that pion cloud effects 
on iV-A transition are significant. 

We present in Fig. 17 our predictions of asymmetry A for several incident electron energies 
E = 0.8, 1.5,2,4 GeV and at the A peak. Here we note that our result for E = 0.8 GeV 
is about 10 % larger than that shown in Fig. 6 of Ref . [38^ . Furthermore, their A decreases 
with Q 2 in striking difference with our predictions. This is perhaps mainly due to the use of 
different N-A form factor, but the differences in treating the non-resonant amplitudes, which 
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interfere with the dominant iV-A transition amplitude, also play some roles. Experimental 
tests of our predictions in Fig. 17 will be useful in examining the extent to which the axial 
N-A form factor Eq.(15), which was determined in our study|3( of (v^, /i~7r) reaction, is 
valid. 

As pointed out in Ref. 0, the form factor Eq.(15) is rather different from the form 
Gn A {Q 2 ) = (1 - ocQ 2 /{f3 + Q 2 ))G A {Q 2 ) with a = 1.21 and (3 = 2 (GeV/c) 2 which was 
determined in earlier works [27]. It is therefore interesting to see how these two form factors 
can be distinguished by the parity violation asymmetry of p(e,e'). This is illustrated in 
Fig. 18. Experiment test of our prediction shown in Fig. 18 will help distinguish these two 
axial iV-A form factors. 

To end this section, We return briefly to the investigation of Quark-Hadron Duality 
presented in section IV. It is also interesting to explore |7| whether the Quark-Hadron Duality 
also exists in the parity violating asymmetry A. This can be done here by comparing the 
predictions of our hadronic model and that of Parton model. The formula for calculating A 
within the Parton Model was derived in Ref.j3^. Keeping only the u and d quarks in the 
considered large x region, one finds for a proton target 

3G F Q 2 2C lu [u(x) + u{x)} - C ld [d{x) + d{x)\ + Y[2C 2u (u(x) - u{x)) - C 2d {d{x) + d{x))\ 



A 



ira2V2 4(«(ar) + u(x)) + (d(x) + d(x)) 

(55) 



where 

Y = i-(i-y) 2 m 
y = |, (57) 

R = — . (58) 

Note that R is the ratio between the longitudinal and transverse total cross sections. We 
use the values of R calculated from the SL model. The variable y depends on incident 
electron energy E and the energy transfer u. The other coefficients in Eq.(55) depend on 
the electroweak coupling constants for leptons and quarks 

1 4 





= 9 e A9v 




= 9l9 d v 




= 9v9 u a 


C 2 d 


= 9v9 d A 



-- + - sin 2 9\ 



2 3 
- - - sin 9 W , 
1 



— - t- 2 sin 9 W , 

i-2sin 2 ^. (59) 

When the radiative corrections within the standard model is included 0, EH C\ u ~ 
-0.1886, C ld ~ 0.3414, C 2u ~ -0.0359, C 2d ~ 0.0265. These values are only slightly dif- 
ferent from those calculated from Eq.(59). 

The asymmetry of p(e, e') predicted from our hadronic model (solid curves) and Parton 
Model (dashed curves) for several typical electron kinematics are compared in Fig. 19. If 
our hadronic model results (solid curves) are confirmed experimentally, the parity violating 
asymmetry in the A region obviously does not show local Quark-Hadron Duality. The 
situation is similar to the results for the spin dependent structure functions shown in Fig. 14. 
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VI. SUMMARY 



The dynamical model developed in Refs. [3, 0, HI (the SL model) has been extended to 
include the weak neutral current contributions for investigating all possible electroweak pion 
production reactions in the region near the A excitation. The main purpose is to examine 
the Quark-Hadron Duality in the neutrino-induced reactions and to explore how the axial 
N-A transition form factor can be determined by using the parity violating asymmetry of 
p(e, e'). The experimental data for testing our predictions can be obtained at JLab and new 
neutrino facilities. 

We have found that the (e, e') structure functions F\ and F 2 predicted by the SL model 
are in good agreement with the recent dataj^, 0| which had verified more quantitatively 
the Quark-Hadron Duality first observed by Bloom and Gilmanj^]. The predicted structure 
functions for (u, e) and (u, v 1 ) processes also show Quark-Hadron Duality to the extent similar 
to what has been observed in (e, e'). Furthermore, we also predict that Quark-Hadron 
Duality should also be seen in all electroweak reactions on the neutron or equivalently 
the isospin 1 = deuteron-like target. These results suggest that the SL model can be 
a candidate hadronic model for developing a theoretical understanding of Quark-Hadron 
Duality within the Standard Model. For example, it will be interesting to explore which 
parts of the predictions from the SL model can be related to the predictions from the Parton 
Model. This is clearly a difficult question to answer and is beyond the scope of this paper. 

We have also investigated the spin dependent structure functions g\ and g 2 of (e, e'). It 
is found that the Quark-Hadron Duality is not seen in the calculated g\ and g 2 , while our 
results for gi and some polarization observables associated with the exclusive p(e, eV) and 
p(e, eV) reactions are in good agreement with the recent data. Experimental data for g 2 in 
the A region are clearly needed. If our predictions are confirmed, we perhaps have more 
complete information for exploring why the Quark-Hadron Duality breaks down in the spin 
dependent structure functions. 

We should also emphasize that our investigations of Quark-Hadron Duality are rather 
qualitative since they are based on the naive Parton Model. To test local Quark-Hadron 
Duality quantitatively, we need to include target mass corrections 2, 47 1 and consider the 
roles of higher twist effects. In particular, the procedures for including these effects in the 
calculations of spin dependent structure functions g\ and g 2 and parity violating asymmetry 
A must also be developed. 

In the investigation of parity violating asymmetry A of the inclusive p(e,e'), we have 
shown that the non-resonant contribution is small at the W = 1.232 GeV A peak and hence 
a precise measurement of A can be used to improve the determination of the axial N-A 
transition form factor. We have also predicted that the parity violating asymmetry A, like 
the spin dependent structure functions g\ and g 2 , does not exhibit Quark-Hadron Duality. 
Predictions for the experiments which can be conducted at JLab have been given. 
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(a) (b) (c) (d) 




FIG. 1: Non-resonant pion production mechanisms of the SL model. See text for explaining how 
these mechanisms are related to the electrmagnetic current, charged current, and neutral current. 
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FIG. 2: The A excitation mechanism. 
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FIG. 3: Structure functions F\ (left) and Fi (right) of (e, e') for the proton target at Q 2 = 1.5 
(GeV/c) 2 . The dashed curves are calculated from using the CTEQ6 parton distribution functions 
at Q 2 = 10 (GeV/c) 2 . The solid curves in the region near £ ~ 0.6 are calculated from the SL 
model. The data are from Liang et al.jloT]. 




FIG. 4: Structure functions F\ (left) and F2 (right) of (e, e') for the proton target. The dashed 
curves are calculated from using the CTEQ6 parton distribution functions at Q 2 = 10 (GeV/c) 2 . 
The solid curves are the results at Q 2 = 0.7, 1.5, 2.5^3.5 (GeV/c) 2 (from left to right) calculated 
from the SL model. The data are from Liang et al.|10j. 
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FIG. 5: Structure functions xF\ (left) and -F2 (right) of (e, e') for the proton (top) and an 1=0 
deuteron-like target (bottom). The dashed curves are calculated from using the CTEQ6 parton 
distribution functions at Q 2 = 10 (GeV/c) 2 . The solid curves are the results at Q 2 = 0.4,1,2,4 
(GeV/c) 2 (from left to right) calculated from the SL model. 





FIG. 6: Structure functions xF\ (left), F2 (center), and xF% (right) of (u, e) for the proton (top) and 
an 1=0 deuteron-like target (bottom). The dashed curves are calculated from using the CTEQ6 
parton distribution functions at Q 2 = 10 (GeV/c) 2 . The solid curves are the results at Q 2 = 
0.4,1,2,4 (GeV/c) 2 (from left to right) calculated from the SL model. 
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FIG. 7: Structure functions xF\ (left), F2 (center), and xFs(right) of (u, u') for the proton (top) and 
an 1=0 deuteron-like target (bottom). The dashed curves are calculated from using the CTEQ6 
parton distribution functions at Q 2 = 10 (GeV/c) 2 . The solid curves are the results at Q 2 = 
0.4,1,2,4 (GeV/c) 2 (from left to right) calculated from the extended SL model described in this 
work. 




FIG. 8: The solid curves are the structure function F2 of p(e, e') as a function of the scaling variable 
£ calculated at Q 2 = 10, 15, 20 (GeV/c) 2 (from left to right) from the SL model. The dashed curve 
is calculated from using the CTEQ6 parton distributions at Q 2 = 20 (GeV/c) 2 . 
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FIG. 9: The beam asymmetry Alt' (defined in Ref.|22|) vs <^>* of p(e, e'ir + )n reaction at Q 2 = 0.4 
(Ge V/c) 2 and W = 1.22 GeV predicted by the SL model are compared with the data of Joo et al. 
[22j . (9* , (/>*) are the pion angles in the center of mass frame of the final ttN system. 
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FIG. 10: The transverse-longitudinal interference term aiT' (defined in Ref . \22\ ) vs cos6* for 
p(e,e'ir°)p (top) and p(e, e'i: + )n (bottom) at W = 1.22GeV predicted by the SL model are com- 
pared with the data of Joo et al. j^] • 0* is the pion scattering angle in the center of mass frame 
of the final ttN system. 
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FIG. 11: The beam-target asymmetries ^(upper half)) and ^4 e f (lower half) (as defined in Ref. 2l|) 
of p(e, e'ir )p reaction as a function of cos 9* integrated over 0* at W = 1.22 GeV and 
0.5 (GeV/c) 2 < Q 2 < 0.9 (GeV/c) 2 (left), 0.9 (GeV/c) 2 < Q 2 < 1.5 (GeV/c) 2 (right). (0*,0*) 
are the pion angles in the center of mass frame of the final irN system. The solid curves are the 
predictions of the SL model. The data are from Biselli et al. [ilj 
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FIG. 12: The beam-target asymmetries ^(upper half)) and ^4 e f (lower half) (as defined in Ref. 2l| ) 
of p(e, e'TT°)p reaction as a function of (jf integrated over cos 0* at W = 1.22 GeV and 
0.5 (GeV/c) 2 < Q 2 < 0.9 (GeV/c) 2 (left), 0.9 (GeV/c) 2 < Q 2 < 1.5 (GeV/c) 2 (right). (0*,0*) 
are the pion angles in the center of mass frame of the final ttN system. The solid curves are the 
predictions of the SL model. The data are from Biselli et al. [ilj • 
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FIG. 13: Spin dependent structure functions gi (left) and 52 (right) for the proton target. The 
dashed curves are from the fits to the deeply inelastic scattering data as explained in the text. 
The solid curves are the results at Q 2 = 0.21,0.35,0.62,0.92,1.37 (GeV/c) 2 (from left to right) 
calculated from the SL model. The data for g\ are from Fatemi et al. [23j. 
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FIG. 14: The parity violating aymmetry terms A ,4 (solid curve) and —Ay (dotted curve) (defined in 
Eqs.(53)-(54)) of p(e, e') as a function of invariant mass W calculated from the extended SL model 
described in this work. The results are for incident electron energy E = 1 GeV and scattering 
angle 6 = 110°. 
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FIG. 15: The parity violating asymmetry terms A^solid curves) and Ay(dotted curves) (defined 
in Eqs.(53)-(54)) of p(e,e') as a function of the incident electron energy E e calculated from the 
extended SL model described in this work. The results are for invariant mass W = 1.232 GeV and 
electron scattering angles 9 = 60° (left) and 9 = 110° (right). 
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FIG. 16: The parity violating asymmetry term (defined in Eq. (54)) of p(e,e') as a function 
of Q 2 evaluated with an invariant mass W = 1.232GeV and electron scattering angle 9 = 110°. 
The solid curve is from the extended SL model described in this work. The dashed curve is from 
turning off the pion cloud effect on the axial N-A transition. The dotted curve near the bottom is 
from keeping only the non-resonant amplitude in the calculation. 
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FIG. 17: The (^-dependence of the scaled parity violating asymmetry A/Q 2 of p(e, e') calculated 
from the extended SL model described in this work. The results are for W = 1.232 GeV and 
incident electron energy E= 0.8 (solid curve), 1.5 (dotted), 2.0 (dashed), 4.0 (dash-dotted) 
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FIG. 18: The scaled parity violating asymmetry A/Q 2 of p(e,e') calculated from the extended 
SL model described in this work. The solid (dashed) curve is from using the axial N-A form 
factor Gff A of Ref-H ( Refs.0). The results are for invariant PF=1.232 GeV and electron angle 
9 = 110°. ' 
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FIG. 19: The scaled parity violating asymmetry A/Q 2 of p(e,e') as a function of the scaling 
variable £ predicted by the extended SL model described in this work (solid curves) is compared 
with that calculated from Eqs.(55)-(59) using the CTEQ6 parton distribution functions (dashed 
curves). The results are for E = 4 GeV ,Q 2 = 0.5 (GeV/c) 2 (top left), E = 6 GeV ,Q 2 = 1 
(GeV/c) 2 (top right), E = 8 GeV, Q 2 = 5 (GeV/c) 2 (bottom left), and E = 12 GeV, Q 2 = 10 
(GeV/c) 2 (bottom right). 
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